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ABSTRACT 

Understanding star formation is problematic as it originates in the large scale dynamics 
of a galaxy but occurs on the small scale of an individual star forming event. This 
paper presents the first numerical simulations to resolve the star formation process on 
sub-parsec scales, whilst also following the dynamics of the interstellar medium (ISM) 
on galactic scales. In these models, the warm low density ISM gas flows into the spiral 
arms where orbit crowding produces the shock formation of dense clouds, held together 
temporarily by their external pressure. Cooling allows the gas to be compressed to 
sufficiently high densities that local regions collapse under their own gravity and form 
stars. The star formation rates follow a Schmidt-Kennicutt Ssfr oc Sg a g type relation 
with the local surface density of gas while following a linear relation with the cold 
and dense gas. Cooling is the primary driver of star formation and the star formation 
rates as it determines the amount of cold gas available for gravitational collapse. The 
star formation rates found in the simulations are offset to higher values relative to the 
extragalactic values, implying a constant reduction, such as from feedback or magnetic 
fields, is likely to be required. Intriguingly, it appears that a spiral or other convergent 
shock and the accompanying thermal instability can explain how star formation is 
triggered, generate the physical conditions of molecular clouds and explain why star 
formation rates are tightly correlated to the gas properties of galaxies. 

Key words: stars: formation - stars: luminosity function, mass function - globular 
clusters and associations: general. 



1 INTRODUCTION 

Star formation is the primary driver of galactic evolution 
and provides the birth place for planetary systems, yet 
it is poorly understood. Observationally, a strong corre- 
lation has long been found between the amoun t of gas 
in the galaxy and the rate at which stars form (|Schmidtl 
1 19591 : lKennicuttlll989l : iKennicutt fc Evans! \20v2i ) . whereby 
the rate of star formation per unit area in the galactic 
disc is correlated with the (total) surface density of gas in 
the disc as Esfr. oc £g a t. This (Schmidt-Kennicutt) rela- 
tion holds both globally in galaxies over a large range of 
galaxy properti es as well as locally o n kpc scales of indi- 
vidual galaxies (|Kennicutt et al]|2007t ). In addition, a near 
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linear relation has been found betw een the dense molecu- 
lar gas and the star formation rat es (Kennicutt et al. Il2007l : 
iBigiel et al.ll200Sl : IWu etafl l201fj| : ISchruba et all 120111 1. al- 
though it has recently been argued that syst ematic varia- 
tions could exist between individual galaxies IShettv et al.l 
(2012). These correlations have inspired many theoretical 
ideas involving everything from gravitational instabilities, 
spiral shocks, local collapse, self-regulation t hrough feedback 
JSill3ll997l; [Elmegre en 2002; Li et al.|[20ul ; iKrumholz et al.l 
l2009l ; IOstriker fc Shettv||201ll ) but until now there has been 
no means of testing these various ideas and thus construct- 
ing a complete theory of star formation in galaxies. 

Recently, deviations from the above star formation re- 
lations have been found on smaller extragalactic (<200 pc) 
scales in the form of large scatter at low gas s urface den- 
sities (|Schruba et al.ll2010l : TOnodera et alj|2010h . These de- 
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viations could be due to intrinsic variation in the cloud to 
cloud star forming properties, sampling issues when rates 
are estimated from integrated fluxes due to massive stars, 
or to the movement of newly formed stars out of their na- 
tal environments. In contrast, star formation rates in nearby 
Galactic molecular clouds, where the rates are calculated di- 
rectly from the young stars present and cloud masse s from 
extinction (E vans et alj|2009 ; iHeiderman et al.ll20ld) . show 
instead a threshold at low surface densities before converg- 
ing to a more sta ndard Schmidt-Kennicutt relation (see also 
lLada et alj|2012h . 

Star formation in galaxies involves the collection of gas 
into dense molecular clouds, and the collapse of localised 
regions of these clouds to form stars, and radiative and ki- 
netic feedback from the young stars back into the interstellar 
medium, processes which may or may not be simultaneous. 
Galactic scale numerical simulations have shown that in spi- 
ral galaxies, molecular clouds likely form thr ough the coa - 
lescence of less dense clouds in spiral shocks i|Dobbsll200sl ). 
aided by self gravity for the most massive clouds. Several 
studies have attempted to study the Schmidt-Kennicutt re- 
lation dr i ven by combinations of s e lf-gra vi ty and feedback 
( LietalJ 120061; IWada fc Normanl 120071; pTasker fc Brvanl 



2009; 



2006; IDobbs fc Pringld 
Dobbs et al.l l2011al; iKim et al. 



Koyama & Ostriker 
201 ll ; lAgertz et all 



2009 



2011 



Shettv fc Ostrikej2012h .' but were unable to resolve the scale 



of star formation and required recipes for star formation and 
sub-grid modelling to account for the combined effects of 
cooling and stellar feedback. 

On the smaller, parsec scales of individual molec- 
ular clouds, local simulations do resolve star formation 



llSmith et al l 120091 ; iBonnell et ail 1201 ll ; iKrumholz et all 
l20ld ; iBatel I2012T) . but these start with highly idealised 
conditions that do not reflect the physics of the cloud 
formation process. Feedback has also been included on 
local scales, but does not appear to significantly alter 
the star formation rates of individual mo lecular clouds 
l|Dale fc Bonnelj|2008l. 1201 ll; [Dale et al.ll2012l). Cloud-cloud 
collis ions l|Whitworth et al.1 Il994l ; iKitsionas fc Whitworthl 
2007) and colliding flows simulation s provide a mecha- 
nism for triggering 



ig star formation jKovama fc Inutsukal 
itj 120061; IVazauez-Sernadeni et al.l 120071; 



2000; Hcitsch et al. 

Baneriee et al.ll2009l ; IClark et al .1120121 ') but lack the galacti c 
context that could generate such flows (|Bonnell et al.tf 2006i. 

This paper presents the first numerical simulations ca- 
pable of following both the large-scale galactic flows of the 
gas and still resolve down to the small scale physics of star 
formation. We use a series of simulations to probe from the 
large scale of a galaxy to the small scale of a star forming 
molecular cloud, in order to investigate what physics may 
lie behind the observed Schmidt-Kennicutt relation. In this 
paper, we concentrate on the physics of cloud formation and 
self-gravity as the possible mechanism behind the galactic 
star formation rates. 
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Figure 2. The equilibrium temperature (broken line) and pres- 
sure is plotted as a function of gas density. The gas is nearly 
isothermal at Ki 10 4 K for low densities, ensuring a high exter- 
nal pressure. Gas above densities of 1 cm -3 cools and is com- 
pressed by the external pressure until reaching densities of greater 
than 100 cm" 3 where pressure balance is restored. Cool gas at 
T < 1000 K can be retained between shocks/collisions as long as 
the external pressure does not drop below P/kg PS 1500 K cm -3 . 



probe three, successively smaller scales, from the scale of a 
spiral galaxy, to the formation of dense clouds in the ISM, 
to the relatively small, sub-pc scale where star formation 
occurs. In addition to following the full galaxy disc simu- 
lation for 360 Myrs, we carried out high-resolution simu- 
lations of a sub-region of the disc with and without self- 
gravity. Our high-resolution non self-gravitating Cloud re- 
simulations followed this region for the last 53.5 Myr of the 
full disc evolution while our Gravity re-simulations including 
self-gravity followed the same region over 5 Myrs near max- 
imum compression to determine the star formation rates. 
Self-gravity is only included in the last set of simulations 
due to the short timestep required to resolve high-density 
self-gravitating gas. Each smaller scale simulation uses the 
conditions from the larger scale to ensure that the global dy- 
namics are included. We neglect feedback from young stars 
and magnetic fields in these calculations in order to probe 
the effects of the spiral shock, cooling and gravity in trigger- 
ing star formation and generating the star formation rates. 
An overview of our fiducial simulations, with a mean surface 
density of 4.2 Mq pc -2 , is shown in Figure 1. In addition, 
we also repeated the Cloud, and then Gravity re-simulations 
at 10 times higher and lower surface densities (discussed in 
§6). All the calculations are detailed in Table 1 and discussed 
further in Section 2.2. 



2 CALCULATIONS 

The numerical simulations used the three-dimensional 
smoothed particle hydrodynamics (SPH) method to follow 
the gas d ynamics in a gala ctic potential including 4 spi- 
ral arms ( Dobbs et alj|2006h . The simulations are used to 



2.1 Thermal Physics 

The thermodynamics of the gas includes compressional and 
shock heating as well as heating from a background inter- 
stellar radiation field. Cooling is included through a cool- 
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Figure 1. An overview of the numerical simulations with E^isc = 4.2Mq pc is shown in terms of their column densities. In the 
global galactic disc simulation, gas shocks, cools and coalesce in the spiral arms to form dense clouds. One of the clouds from the global 
simulation is used to perform the Cloud re-simulation, where we have increased the mass resolution by a factor of 256 in order to study 
in detail how the star formation region is formed. The last 3 panels show the cold gas (T < 50 K) in the central regions of the Gravity 
re-simulation, which now includes self gravity and sink particles to indicate stellar clusters. Gas in the two re-simulations explores scales, 
and exhibits the densities and temperatures typical of molecular clouds (10 3 cm -3 and < 20 K) which are otherwise unfeasible in full 
global disc simulations. In each panel, the column density of the gas is shown, the logarithmic scale from 10~ 4 to 0.1 g cm -2 for the 
global simulation, 10 — 3 to 0.1 g cm -2 for the Cloud re-simulation, and from 0.002 to 2 g cm -2 for the Gravity re-simulation. 
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Model 


^disc 


Particle mass 


Region size 


Start time 


Total time 




(M Q pc- 2 ) 


(M ) 


(PC) 


(Myr) 


(Myr) 


Global simulation 


4.2 


40 


20,000 





370 


Cloud re-simulation 


4.2 


0.156 


2000 


310 


54 


Gravity re-simulation 


4.2 


0.156 


250 


350 


5 


Cloud re-sim high 


42.4 


1.56 


2000 


345 


17 


Gravity re-sim high 


42.4 


1.56 


250 


350 


2 


Cloud re-sim low 


0.42 


0.0156 


2000 


345 


10 


Gravity re-sim low 


0.42 


0.0156 


250 


350 


8 



Table 1. This table displays the details of the simulations presented. The times shown are in relation to the global simulation. The 
Gravity and low and high re-simulations are all derived from the Cloud re-simulation. The total masses are 2.5 XlO 9 Mq for the global 
disc simulation, 1.71 XlO 6 for the Cloud re-simulation and 1.61 xl0 B for the Gravity re-simulation. For the low and high re-simulations 
the total masses are ten times lower (fa 10 5 Mq) and higher (sj 10 7 Mq) than for the Cloud and Gravity re-simulations, respectively. 



ing function appropriate for interstellar gas that is optically 
thin, including atomic and molecular lin e cooling and cool- 
ing from dust (|Kovama fe Inutsukal2002l ) . This cooling func- 
tion is used to calculate an equilibrium temperature and a 
cooling time to reach this equilibrium. Together this allows 
the internal energy to be integrated implicitly and without 
aff ecting the hydrodynamical tim esteps (for more details, 
see IVazquez-Semadeni et al.l l2007). At low densities, n < 1 
cm -3 , the gas is nearly isothermal at T ~ 10 4 K. At higher 
densities, a cooling instability occurs when atomic line cool- 
ing dominates at densities of n > 1cm -3 while molecular 
and dust cooling is more important a t the high gas densi- 
ties a ppropriate for molecular clouds l|Kovama fe Inutsukal 
2002). The cooling curve has equilibrium temperatures that 
produce a bi-stable ISM where a warm (T « 10 4 K) low 
density medium exists in pressure equilibrium with the 
cold (T k, 1QK) dense gas where sta r formation can occur 
i|Wolfire et al.ll 19951 ; [Field et alj|l969h . This parametrisation 
neglects the detailed chemistry and locally variable extinc- 
tion and will likely overestimate the cooling rates at low den- 
sities (and overestimating the generation of self-gravitating 
gas and the star format ion rates) while underestimating 
them at higher densities (|Clark et al.ll2012l ; iGlover fe Clark! 
l2012al ). 

The critical aspect of the simulations is the thermal 
instability that occurs when the gas is compressed to den- 
sities > 1 particle cm -3 , and is best understood in terms 
of the gas pressure. Figure [5] shows the equilibrium pres- 
sure density relation for this thermal cooling law. The gas 
enters the shock at low densities and is compressed nearly 
isothermally at w 10 4 K, increasing its internal gas pressure. 
When gas is compressed beyond a critical density of 1 cm -3 , 
cooling dominates and the gas cools and is under-pressured 
until reaching densities above 100 cm -3 , the densities of 
star forming molecular clouds. The surrounding gas ensures 
a high external pressure and along with the ram pressure 
of the shock, compresses the cooling gas to densities where 
self-gravity can take over. 

Once the gas leaves a shock, the external pressure is 
reduced and the gas can return to the warm phase. If a 
sufficiently high pressure is maintained, then some of the 
gas can be retained in a cool phase at densities near n « 10 
cm" 3 and temperatures of 100 < T < 1000K. This is likely 
to occur between collisions while the gas is still in a given 



spiral arm. Subsequent collisions involving the cool gas will 
be stronger due to the now lower sound speed in the gas. 



2.2 Numerical Simulations 

The global disc simulation used 2.5 x 10 7 SPH particles to 
model 10 9 Mq of gas in an annulus 5 to 10 kpc, for an 
initial surface density of Edisc ~ 4.2 Mq pc -2 . The gas was 
given an initial temperature of 10 4 K, and allowed to cool 
during the evolution due to the thermal physics, resulting 
in a multiphase ISM. The evolution is followed for over 350 
Myrs to provide the initial conditions and boundary regions 
for the high resolution simulations. 

The high resolution re-simulations were constructed by 
taking a high-density region from the global disc calculation, 
at a time of 350 Myr (3.5 x 10 8 years), and re-resolving this 
by splitting the original particles into 256 lower-mass par- 
ticles. A 250 pc region of interest was chosen by eye, con- 
taining 1.71 x 10 6 Mq of cold gas, and then modeled using 
1.09 xlO 7 SPH particles with an individual particle mass of 
0.156 Mq. The region was traced backwards through ~ 1/4 
of a galactic orbit (~ 5 x 10 7 years) to produce the initial 
conditions for this high resolution re- simulation. The region 
was augmented to include boundary particles, regular SPH 
particles but outwith our region of interest, to ensure the hy- 
drodynamics of our re-simulation are accurately modelled. 
These boundary particles include any particle that comes 
within 40 % of the size of the region. Any of these boundary 
particles that come within 15 % of the size of the region 
are split into 4 lower-mass particles, while those that come 
within 10 SPH kernel lengths of a particle in our region 
of interest are split into 16 lower-mass particles. Further- 
more, if these particles come within 4 SPH kernel lengths 
(defined from the original simulation) , they are split into 64 
lower-mass particles. This ensures that the particles in our 
region of interest, and all other particles in the simulation, 
only encounter other particles that are at most a factor of 
4 different in mass. The boundary particles are evolved as 
are the standard SPH particles. The rest of the galactic disc 
was not modeled in these Cloud, and Gravity re-simulations. 
The Cloud re-simulation (region of interest plus boundaries) 
contained and 2.8 x 10 7 SPH particles and 3.09 x 10 7 M . 
The simulation was then evolved forwards for 54 Myr as the 
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Figure 3. The evolution of the high-resolution, non self-gravitating Cloud simulation is shown in panel a at 4 different times ( 307.7, 
325.6, 333.1 and 353.8 Myr), as the gas rotates clockwise around the galaxy. The column densities are plotted logarithmically from 10 — 4 
to 0.1 g cm— 2 (bottom scale) and the panel shows a region 12 kpc on a side. The smaller panels (b-f) show the zooms of the evolution, 
from top left to bottom right, at times b) 17.9 Myr, c) 25.4 Myr, d) 33.0 Myr, d) 44.1 Myr and e) 53.5 Myr from the start of the 
high-resolution simulation (307.7 Myr). Each panel shows a 200 by 200 pc region centred on the clump as it forms though the spiral 
shock and cloud coagulation. The gas is compressed into individual clouds which interact with each other and the inflowing cold gas. 
After maximum compression, the cloud is torn apart by the divergent flow as it leaves the spiral arm. The column densities are shown 
logarithmically between a maximum of 0.1 and a minimum of 0.001 g cm" 2 (top scale). 



gas flowed back through the spiral shock to reproduce the 
region of interest at high resolution. 

The Gravity re-simulations started using the inner 250 
pc of the Cloud re-simulations at a ti me of 350.8 Myr. Self- 
gravity is calculated via a tree-code (|Benz et alj|l990h and 
includes both the galactic potential and the gas contribu- 
tions. The particle mass of 0.156 Mq in the Gravity re- 
simulation is adequate to resolve gravitational fragmenta- 



tion down to masses of ~ 11 Mq. This is not sufficient to 
resolve the formation of individual low mass stars but is cer- 
tainly sufficient to resolve the formation of stellar clusters 
where most stars form. These values are ten times lower and 
higher in the low and high re-simulations, respectively. 

Star formation is modelled by sink particles (Bate ct al. 

Il995h where bound, collapsing regions are replaced by a sin- 
gle particle that can continue to accrete infalling gas that 
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falls within 0.25 pc. Replacing a full SPH kernel containing 
some 70 SPH particles ensures that the star formation is 
resolved and driven by the local gravitational collapse. This 
allows us to derive the star formation rates rather than input 
a parametrised rate based on local properties as is commonly 
used in cosmological simulations. We assume a 100 % effi- 
ciency of star formation for the gas that falls within our sink 
particles. Star formation efficiencies within 0.25pc could be 
as low as 10-50 %. and hence we will determine an upper 
limit on star formation rates by a factor of 2-10. 

For the two additional Cloud and Gravity re- 
simulations, the mean surface densities in the disc were in- 
creased and decreased by an order of magnitude (Edisc = 
0.42 and 42 Mq pc -2 ) compared to the standard Edisc = 4.2 
Mq pc -2 run. The global dynamics have minimal depen- 
dence on the different gas densities in the absence of self- 
gravity. However the gas is slightly warmer/cooler in the 
low/high surface density re-simulations due to the depen- 
dence of heating and cooling on density. We thus reran the 
non self-gravitating Cloud re-simulations for a minimum of 
10 Myrs to ensure the gas temperatures were correct. The 
Gravity re-simulations were evolved for shorter and longer 
time periods for the higher and lower surface densities, re- 
spectively, reflecting the different free-fall timescales of the 
self-gravitating gas involved. In all cases, the star forma- 
tion rates had reached a steady level before the end of the 
simulation. 




i_d i i i i i i i i I 
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Figure 4. Column densities of the different g function 
of temperature are shown at maximum compression in the non 
sclf-gravitating Cloud re-simulation. The hot (T ^ 5000K) and 
warm (500 ^ T ^ 5000 K) gas is plotted between minimum and 
maximum column densities of 1 X 10 -4 and 5 X 10 -3 g cm -3 
(top scale). The cool (50 s£ T sC 500K) and cold (T sg 50K) gas 
is plotted between minimum and maximum column densities of 
2 X 10~ 3 and 1 X 10 _1 g cm 3 (bottom scale). 



3 RESULTS 

3.1 Global disc simulation 

The global disc simulation follows the gas, with an initial 
surface density of Edi sc ~ 4.2 Mq pc -2 , through 350 Myr in 
which time it undergoes several (3-4) spiral arm passages. 
During the simulation, the gas is subject to the spiral po- 
tential, heating and cooling, but not self-gravity of the gas. 
The spiral arm potential produces con vergent gas stream s 
that meet and shock in the spiral arms (|Dobbs et alj|2006l) . 
Each time the gas enters the spiral arms, the gas is com- 
pressed, cools and then re-expands in the inter-arm region. 
The cooling instabilities in the spiral arms lead to clumpy 
molecular cloud structures, which as the gas leaves the arm, 
are sheared into spurs or feathers In addition, there are re- 
peated collisions and internal shocks inside the spiral arms. 
This ensures that the structures we are studying are formed 
self-consistently, and thus provide an appropriate starting 
point for studying how star formation is initiated in galactic 
discs. 

At the end of the global disc simulation (Fig. 1), pa 15 
% of the gas mass (1.5 X 10 8 M Q ) is sufficiently cold (T < 
100K) and dense (n > 100 cm -3 ) to provide the necessary 
environment for star formation. The thermal instability and 
the repeated compression generates a multiphase interstellar 
medium where cold, dense g as is in pressure e quilibrium with 
the warm, low density gas (|Field et al.lll96i^ ). 



3.2 Cloud re-simulation: Formation of Molecular 
Clouds 

Figure [3] shows our standard (E ~ 4.2 Mq pc -2 ) high res- 
olution Cloud re-simulation as the gas evolves around the 
galactic disc over 54 Myr, along with several zoomed snap- 
shots of the gas in its centre of mass frame. This shows 
how the gas clouds form through infall in the spiral shock 
and cloud coagulation in the spiral arm. The cloud complex 
consists of material which has been in the spiral arm for 
over 10 Myr and gas which has entered from the inter-arm 
region. Gas from the spiral arm is already cold and dense, 
whilst gas entering comprises warm (> 8000K) gas and some 
cooler and denser gas from previous spiral arm encounters. 
These large scale flows of gas converge to produce a complex 
structure of dense gas clouds including a central (~ 40 x 80 
pc) cloud (see Figure 1) . The clouds are not well defined ob- 
jects as material is continually entering and leaving them. 
The cloud is eventually shredded by the galactic tidal field 
after 53.6 Myrs. 

At each shock or collision, the gas is compressed and 
cools while maintaining pressure equilibrium with the sur- 
rounding warm gas. The cool gas is embedded in the warm, 
low density phase (Figure [4| and subsequent shocks involve 
this clumpy medium. If the external pressure decreases (i.e. 
after a shock), the cool gas can re-expand. If the pressure 
drops below the critical level to maintain a bistable gas, 
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Figure 5. The relation between the mean 1-D velocity dispersion 
and the sizescale of subregions in the cloud centred on the dense 
gas is shown from the early stages of the Gravity re-simulation. 
The rms scatter in the data is indicated by the errorbars. Points 
without errorbars represent only single global values are possible 
at that size. The relation follows the typical Larson relation for 
turbulent giant molecular clouds from 0.1 to 10 pc, but with sig- 
nificant scatter. The Larson relation Sv = 0.9Rp C 5 (km/s) is also 
plotted as the solid line. 



Figure 6. The relation between the binding mass needed to be 
gravitationally bound and the enclosed gas mass is plotted for 
various regions centred on local high densities throughout the 
Gravity re-simulation. The initial values are plotted as red tri- 
angles. Values determined for the same cloud but at increasing 
radii are plotted individually. A ratio of 1 (bound) and 10 (highly 
unbound) are plotted as solid lines. Only on small scales, and 
hence low masses, are the clouds bound even when self- gravity is 
included. 



all the cool gas returns to the warm phase. If the external 
pressure is instead maintained above this critical level, then 
some of the gas will remain cool. This is especially true for 
the multiple collisions inside a spiral arm. 

In this Cloud re-simulation, approximately half of the 
mass (excluding boundary particles) reaches molecular cloud 
densities and temperatures. The high density regions are 
surrounded by a cocoon of warm gas which provides the 
mass reservoir from which they form, but crucially also helps 
hold them together. The cold dense gas occupies a hole in 
the warm gas (Figure [4]). 

The high velocity shock and cooling also produce 
a turbulent and highly structured env ir onment (see also 
Bonnell et all l200d : IHeitsch et all 120061 : iDobbs fc Bonnelll 
20071 ). The turbulence provides an important global support 
to the clouds, while the structure provides the seeds for local 
gravitational collapse. The dense gas leaves the spiral arm 
some 6 Myr after maximum compression and is torn apart 
by the galactic shear , setting an upper limit on the lifetimes 
of molecular clouds (|Dobbs et al.ll2006T ). 



3.3 Gravity re-simulation: Following Star 
Formation 

In the Gravity re-simulations, the inner 250 pc radius re- 
gion of the Cloud re-simulations were reran near the max- 
imum compression including self-gravity to follow the star 
formation process. Figure 1 also shows three frames of the 
central regions from our standard Edisc ~ 4.2 Mq pc -2 
Gravity re-simulation, as seen from inside the plane of the 



galaxy. The clouds are very structured with local star forma- 
tion occurring along the filaments while large clusters form 
where the filaments intersect, and are similar to star form- 
ing regions observed by the Herschel in frared space telescope 
jMolinari et al.ll2010l ; iHill et all 1201 ll ). The clouds are not 
well defined objects with clear boundaries. Instead, the cold 
dense gas represents the peaks in a continuous distribution. 
Furthermore, gas enters and leaves these 'clouds' throughout 
their lifetimes. 

Star formation commences approximately 0.1 Myr af- 
ter the start of the Gravity re-simulation, and continues 
for the full 5 Myr evolution. Local regions that become 
self-gravitating have sizes of a few pc and thus densities 
> 100 — 1000 Mq pc -3 . These regions are colder in their 
centres with minimum temperatures reaching < 20K. The 
total star formation rate reaches a few times 10~ 2 Mq yr _1 . 
The star formation efficiency, the fraction of the gas turned 
into stars, is a few per cent in the first few million years and, 
in the absence of feedback and magnetic fields, increases to 
over ten per cent after five million years. 

In addition to producing the filamentary and clumpy 
structure in the gas, the shock drives large scale 'turbu- 
lent' motions into the dense gas. The incoming clumpy 
gas combined with the convergent flows in the spi- 
ral arm generate a large velocity dispersion ( see also 



Bonnell et al ]|2006l; IHeitsch et al ] |2006l ; lKovama fc Inutsukal 



20021 ; IDobbs fc Bonnelll 120071 ). Figure [5] plots the velocity 
dispersion sizescale relation in the dense gas from near the 
beginning of the Gravity re-simulation. The velocity disper- 
sion is calculated as the mean of many different subregions 
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of a given size centred on the densest gas particles. Indvid- 
ual subregions are calculated at increasing sizes until they 
meet another subregion or an edge in the mass distribution is 
reached (when M oc R n with n < 0.75). These motions fol- 
low the standard Larson-relations where the velocity disper- 
sion sca l es as the square root of the size-scale of the re gion 
jLarsonl 1 198 it iHever fc Brum] |2004| ; iHever et al l 120091 ). In 
fact the simulation reproduces both the approximate slope 
and normalisation (Sv = 0.9i?p C 5 km/s) in the cold dense gas 
up to ~ 10 pc. Figure [5] also shows the rms scatter in the 
data, plotted as errorbars showing that there is a significant 
variation in the local velocity dispersion. 

The velocity dispersion occurs as clumps entering 
the shock on semi-random trajectories encounter differ- 
ing a mounts of gas and hen ce various levels of decelera- 
tion jDobbs fc Bonnellll2007l ). Thermal instabilities in the 
shocked gas also aid in generating the observed velocit y 
dispersions (|Heitsch et al.|[2006l ; iKovama fc Inuts uka 2002). 
The velocity dispersion relation from the simulations be- 
comes flat above ~ 10 pc, in contrast to that observed in 
GMCs. This could indicate that supernova or other feed- 
back, or alternatively a magnetic coupling from the larger 
scale flows, are needed to energise the ISM on larger scales. 

The kinematic motions induced in the shock limit the 
onset of gravitational instabilities in the dense gas. Figure [6] 
shows the relation between the cloud mass and the gravi- 
tational binding mass, the mass needed to be gravitation- 
ally bound given the internal kinematics. We evaluate the 
cloud properties by integrating outwards from the multiple 
dense regions and evaluate the total mass, kinetic energy 
and hence binding masses at each radius. The kinetic ener- 
gies are calculated in the centre of mass frame at each scale 
so that bulk motions are excluded. The resultant properties 
show that the clouds are bound only on small sizescales and 
masses (masses below 100 to 1000 M©) and it is at these 
scales that we expect gravitational instabilities to induce 
star formation. At large radii, the clouds tend to be unbound 
with binding masses of order a few to more than 10 times 
the cloud masses. This is broadly consistent with ob serva- 
tions of Galactic molecular clouds l|Hever et al.ll2009l ). This 
excess kinetic energy at large scales can help maintain low 
star formation efficiencies as significant p arts of the clouds 
may never become gravitationally boun d (jClark et al . 2008; 
iBonnell et aHl201ll : [Dobbs et al.ll2011bh . The external pres- 
sure, including the ram pressure of the incoming gas which 
is responsible for forming the clouds, is an important con- 
tributor on large scales. 

Although the clouds appear to be gravitationally un- 
bound, there is significant collapse in the vertical direction 
throughout the Gravity simulation. The turbulence induced 
by the shock is anisotropic and is significantly lower in the 
vertical direction than in the plane of the galaxy. If sup- 
ported only by the vertical motions, the dense gas would be 
gravitationally bound at all scales. Feedback, magnetic fields 
and any hydrodynamical instabilities that may be missed by 
SPH could act to increase the vertical turbulence and hence 
decrease the star formation rates and efficiencies further. 



4 THE SCHMIDT-KENNICUTT RELATION 

One of the strongest observed relations of star forming galax- 
ies is the Kennicut-Schmidt relation, where the surface den- 
sity of star formation is c orrelated with the surface den- 
sity of gas as Ssfr oc S^ a t l|Schmidtll 19591 ; lKennicuttlll989l ; 
iKennicutt fc Ev ans 2012). In Figure we plot the local rate 
of star formation (SFR), plotted in cells of 50 by 50 pc, 
against the total surface density of the gas in the Gravity 
re-simulation, as seen from within the plane of the galaxy. 
The star formation rates are calculated directly from the 
mass incorporated into our sink particles over the timescale 
of the simulation. 

The star formation rates follow a Schmidt-Kennicutt 
type Espr oc S* a 1: relation whe n considering the tota l 
gas surface density (|Schmidtj|l959l ; IKennicutt etall 120071 ) . 
Included in Figure [7] are the results for two other self- 
gravitating re- simulations, with initial surface densities a 
factor of ten higher (42 M Q pc" 2 ) and lower (0.42 M pc" 2 ) 
than the initial run, respectively. Together there is a strong 
continuous Esfr oc E" as type relation across all three Grav- 
ity re-simulations, arising solely due to the combination of 
the spiral shock, associated thermodynamics and the self- 
gravity of localised regions. A best x-squared fit to the data 
finds n « 1.55 ± 0.06 when all non-zero values of the star 
formation rate are considered over the three Gravity simu- 
lations. If we exclude the low surface density simulation, or 
fit only data with E gas ^ 10 Mq pc" 2 , where a star forma- 
tion threshold may be present, we find a best fit slopes of 
n = 1.42 ± 0.07, and n = 1.44 ± 0.05, respectively. 

In contrast to the success in obtaining a reasonable 
slope for our star formation relation, the normalisation of 
the data is significantly higher than typical observed values. 
The star formation rates are offset to higher values from 
the Schmidt-Kennicutt relation, althoug h they are closer to 
those found in nearby molecular clouds l|Evans et al.ll2009l ; 
iHeiderman et af1l2010l ) which also use direct estimates of the 
star formation rates. Averaging over larger regions could re- 
duce the relation if it included significant areas with gas 
but no star formation. This is not the case in our present 
simulation as it was deigned to resolve the high-density star 
forming region. It is likely that additional physics is required 
to globally reduce the star formation rates. 

Many surveys have shown that the star formation law 
is shallower when looking at m olecular gas with a di- 
rect Ssfr oc E 

molcc ular 

relation (|Bigiel et al.l 120081 . 1201 it 
ISchruba et al.ll201ll ). This suggests that the star formation 
rates are linked to the formation of dense molecular gas and 
that thereafter there is a direct correspondence at a con- 
stant rate given by the free- fall time and an efficiency factor 
Esfr = eEmoiccuiar/te (e.g. iKrumholz et aT1l2009l ). We also 
find (see Figure[S]a near-linear relation of the star formation 
rates to the cold (T < 100K) and dense (p > 10M Q pc" 3 ) 
gas Edonso- This occurs as the free-fall time of the cold and 
dense gas is roughly constant at a few to 10 Myrs. We expect 
that feedback and magnetic fields or an increase in the ver- 
tical component of turbulence should increase this timescale 
for star formation and decrease the star formation rates in 
a global fashion. 

If the star formation rates follow a linear relation with 
the dense, cold, molecular gas, then it is at the phase where 
molecular clouds are formed that the non-linear nature of 
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Figure 7. The local star formation rates per unit area (left) and the surface density of cold dense gas (right) are plotted against the 
surface density of gas for the three Gravity re-simulations with initial surface densities of 0.42 Mq pc — 2 (blue stars), 4.2Mq pc — 2 (black 
filled triangles), and 42.0Mq pc — 2 (red squares). The points represent local 50x50pc regions of the disc as viewed from the plane of the 
galaxy. The SFRs (left) correspond to a time 3.5 X 10 6 years from the start of the Gravity re-simulations, except for the high-surface 
density run which corresponds to 1.3 X 10 6 years. The cold, T < WOK, and dense, p > lOM© pc — 3 , gas (right) is m easured at the start 
of the Gravity re-simulations. The dashed line is the approximate location of the Schmidt-Kennicutt SgpR oc Sgat l|Kennicutt &c Evand 
2012) relation (left) and shows a Edense ok Sg a * relation (right). The dotted and solid lines at the left show two x~squared fits to the data 
representing the total dataset (slope = 1.55 ± 0.06) and for just the two higher surface density simulations (slope = 1.42 ± 0.07). Null 
values are set to 1 0~ 4 Mp yr —1 kpc~ 2 and M pc — 2 , respectively. The small cyan and green (small open triangles) points represent the 
iBigiel et al.ll201(j| and iHeiderman et al.ll2010l Table 1 data, respectively. Note that the Edense oc relation implies a molecular star 

formation relation of Ssfr <x Edense- 



the Schmidt-Kennicutt relation must arise. In Figure [7] we 
also plot the surface density of cold (T < 100K) and dense 
(p > IOM0 pc -3 ) gas Edense, that corresponds to molecu- 
lar clouds, as a function of the total surface density of gas 
as viewed horizontally in the plane of the galaxy. We find 
a Edcnsc oc Eg a t relation with a similar slope as the star 
formation Schmidt-Kennicutt relation. This shows that the 
non-linear slope in our simulated S-K relation arises in the 
formation of the cold gas. The 'molecular' gas fraction is 
zero below surface densities of a few Mqpc -2 and increases 
linearly with the surface density from a few percent to being 
pre dominantly cold and d ense at high surface densities (see 
also IKrumholz et al.l l2009h . 

Unexpectedly, the Edense oc E ga | relation is found in 
both the self-gravitating and the non self-gravitating simu- 
lations. That this fundamental relation arises independently 
of self-gravity shows that gravity cannot be the determining 
physics in generating the cold gas, and ultimately in setting 
the star formation rates. Instead, it is the shock and thermal 
physics of how the cold dense gas is genera ted which drives 
the re sul ting star formation ra t es (see also Idover fc Clarkl 
l20T2bl lal; IKrumholz et al.ll201ll ; iKrumholj [20121 ). Gravity's 
role is then limited to the collapse of the dense gas (which 
in our simulations has a near constant free-fall time, iff) such 
that Esfr ~ Edense /iff oc Edense- A slower collapse due to 
magnetic fields or equivalently a relaxation of our assump- 
tion of 100 per cent star formation efficiency would then 
explain the offset found in the star formation rates. 



The steeper relation between the dense gas and the to- 
tal gas present must arise from the shocks and accompany- 
ing cooling. In the re-simulations it is found that the cold 
dense gas is formed by collisions between clumps of cool 
(100 < T < 1000 K) gas formed in previous shocks. Fig- 
ure shows the evolution of gas that comprises a region of 
high star formation rate over the 50 Myr of the Cloud re- 
simulation. The gas spends only a relatively short time (5 
Myr) as dense and cold (T < WOK ; p > 10M Q pc" 3 ) but 
is mostly cool (T < 1000 K) throughout the full 54 Myr of 
the re-simulation. This 'cool' gas (often termed the warm 
neutral medium WNM) is due to earlier shocks during the 
spiral arm passages and for ms a clumpy mediu m embedded 
in the warmer gas (see also lDobbs et al]|2012h . Subsequent 
collisions involving this cool gas result in much stronger com- 
pressions, hence the formation of the cold, dense (T < WOK ; 
p > IOM0 pc -3 ) gas and ensuing star formation. Where the 
probability of an individual clump colliding with another is 
low, this results in a non-linear scaling between the incoming 
E gas and shock- produced Edense- We construct below a sim- 
ple clumpy shock model based on the requirement of clumps 
to interact in order to form the dense gas. This toy model 
can reasonably reproduce the Edense oc E ga t relation when 
clumps are in pressure equilibrium with the gas and their 
sizes are then given by J? oc Eg a i . 
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Figure 8. The local star formation rates per unit area are plot- 
ted against the surface density of cold (T < 100) and dense 
(p > IOMq pc — 3 gas for the three Gravity re-simulations with 
initial surface densities of 0.42 Mq pc -2 (blue stars), 4.2Mq pc -2 
(black filled triangles), and 42.0Mq pc~ 2 (red squares). The 
points represent local 50x50pc regions of the disc as viewed from 
the plane of the galaxy. The SFRs correspond to a time 3.5 X 10 6 
years from the start of the Gravity re-simulations, except for the 
high-surface density run which corresponds to 1.3 X 10 6 years. The 
star formation rates follow a quasi-linear relation with the surface 
density of dense gas (corresponding to molecular gas) consistent 
with EgpR = SmoiecularAff an d a free-fall time of 1 to 10 mil- 
lion years. These simulations do not include feedback or magnetic 
fields which would slow the star formation process. 



4.1 Star formation rates: thresholds 

Both the star formation rates and the surface density of cold 
gas (Fig. [7]) exhibit minimum surface densities below which 
there is no dense gas or star formation. These thresholds are 
due to a minimum external pressure necessary to maintain 
the cool (WNM) gas, from which both the dense gas and 
the star formation activity are derived. At surface densities 
lower than E gaa « IMq pc -2 , the dense gas is over-pressured 
relative to the external medium and expands. It then heats 
up as the decreasing cooling rates at lower densities no 
longer balance the external radiation field. This threshold 
in maintaining the cool gas produces a steep Esfr — E gas 
over a short range. The higher threshold in the 42Mq pc~ 2 
resimulation reflects simply a scaled up minimum surface 
density where cool gas is present in the initial conditions. 

The thresholds occur at higher surface densities in the 
re-simulation that started with a higher mean surface den- 
sity. This is likely to be due to the resolution and the 
methodology used to construct the re-simulations. The ini- 
tial amount of cool gas present in each re-simulation was 
taken from the global disc model and then scaled accord- 
ingly. For the lower surface density re-simulation, this will 
have produced an overabundance of cool gas but this will 
quickly dissipate as it will no longer be in pressure equilib- 
rium with the lower external pressure. In contrast, the higher 



Figure 9. The mass fraction of cool (T < 1000 K, dashed line), 
cold (T < 100 K, solid line) and dense gas (p > 10M Q pc -3 ) 
gas is plotted for a region of high star formation throughout the 
Cloud re-simulation. The gas spends only a short time ( 5 Myr) 
as cold and dense and hence capable of undergoing gravitational 
collapse and forming stars. In contrast, the gas is mostly cool 
(T < 1000 K) throughout the evolution before the spiral shock 
compresses the cool gas t o form the dense cold gas where star 
formation occurs (see also lDobbs et alj|2012b . 

surface density re-simulation will start with a deficit of cool 
gas. Given a sufficiently long evolution, this cool gas would 
be self-consistently regenerated but it is apparent that over 
the 50 Myrs of the re-simulations, this has not occurred. 
Thus, there is a physical threshold that emerges from the 
simulations but in the high surface density re-simulation has 
been scaled to higher surface densities due to the resolution. 



5 STAR FORMATION RATES: A CLUMPY 
SHOCK MODEL 

The star formation rates scale directly with the surface den- 
sity of cold gas. We therefore need to understand the re- 
lation between the dense gas and the total surface density 
Sdense oc E ga t . As this occurs even in the absence of self- 
gravity, we require a model that includes only the shock 
compression and cooling. In our simulations the cold gas is 
primarily produced by the collision of cool clumps of gas 
which are compressed to higher density cold gas by the ad- 
ditional ram pressure. The likelihood of cold gas formation is 
therefore determined by the probability of cool gas collision. 

The amount of cold, dense gas generated in a shock 
Sdense * s gi ven °y the amount of cool gas entering the shock 
and the probability of an individual clump of cool gas hitting 
another clump of cool gas, so 

Edcnsc OC EdP (1) 

where E c i is the surface density of cool clumps and P is 
the probability of each clump colliding with another. Let us 
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Figure 10. The surface density of cool gas (T < 1000 K) rela- 
tive to the total gas density for the Cloud Re-simulations with 
mean surface densities 0.42 Mq pc — 2 (blue stars), 4.2Mq pc -2 
(black filled triangles), and 42.OM0 pc — 2 (red squares). Above 
the threshold needed to sustain cool gas the surface density of 
cool gas scales almost linearly with the total surface density. 



consider the maximal and minimal cases which the true col- 
lision probability, P, must lie between. In the regime where 
the clumps occupy a large fraction of the total volume, all 
clumps will hit another one, and the rate of forming 

dense gas from these collisions will scale as E c i. In the regime 
where the clump volume fraction is low, the probability of 
a given clump colliding with another is then E c i7ri? 2 , where 
-kB? is the cross section for collision of a clump of radius R. 
If R is independent of environment, then the rate of form- 
ing dense gas from the clumpy shock will scale as E 2 !. As 
the true clump collision probability lies between these two 
regimes, the formation of dense gas must scale as the cool 
gas surface density to a non-linear power between 1 and 2. 

To estimate the value of the true collisional rate we 
consider the following simplified toy model. Firstly we as- 
sume that the cool gas clumps are in pressure equilibrium 
with their environment, and that all the clumps have a 
similar gas temperature. This means that the radius of a 
given clump, R, will be proportional to its surrounding pres- 
sure, p, which is proportional to the gas surface density as 
R oc p -1 / 3 oc T,gal 3 ■ Putting this into our above relations 
this gives 



oc E„iE 



-2/3 



cl 'c/as 



(2) 



Figure [TOl shows the surface density of cool clumpy gas 
verses the total gas density. As previously discussed, below 
a threshold density of around 1 Mqpc -2 there is insufficient 
pressure to confine the cool gas, however above this thresh- 
old the surface density of cool clumpy gas scales close to 
linearly with the gas density, E c ; oc E gaa . The dense gas 
surface density then becomes roughly 




Figure 11. To test the probability of cold gas generation from 
collisions we randomly place circles of radius R oc Sg a ™ in a given 
area and calculate the surface density of dense gas where they 
overlap. 




oc E 



4/3 



(3) 



100 1000 io 4 



Figure 12. The resultant surface densities (in dimensionless 
units) of dense (shocked) gas versus total gas for the simple 
clumpy shock model where R oc E" 1 / 3 . A clumpy shock model 
can reproduce a reasonable relation between the surface density 
of dense gas and total surface density provide R oc S — n with n 
in the range of 1/3 and 0.4. The dashed line shows a slope of 1.4. 



in close agreement with the relationship found in our full 
simulation. 

To further investigate these ideas we make the following 
test of the probability of randomly generated clumps over- 
lapping. In a given area we randomly put down clumps with 
a radius R oc E~ a ™. Where clumps overlap we assume cold 
dense gas will be produced and compare its surface density 
with the total cool gas surface density. Figure [11] shows a 
schematic of this process. We repeat the test varying n and 
find reasonable Edense oc Eg a 4 s relations when n is between 
1/3 and 0.4. Figure IT2"1 shows the relationship for n = 1/3. 



6 CONCLUSIONS 

The results presented here show that a spiral shock coupled 
with the cooling rates in the interstellar medium can ex- 
plain the origin of the Schmidt-Kennicutt relation for star 
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formation rates in galaxies. The shock formation of cold gas 
provides the link between the total gas surface density and 
the ensuing star formation rates. It also explains the linear 
relation found between the dense molecular gas and the star 
formation rates, such that E 

SFR 2_/dense C< Z-/o- as . 

Self-gravity is only required to operate locally in the 
dense gas which has a characteristic free-fall time. The pres- 
ence of cooler gas (T < 1000K) in pressure equilibrium with 
the warm gas as it enters the shock is important in forming 
the den se, cold self-gravita ting clouds where star formation 
occurs l|Pringle et alj|200lh . A minimum pressure or surface 
density is required to maintain the cool gas, and hence sets 
a threshold for star formation. A simple shock model for the 
pressure-bounded cool gas can account for the non-linear 
nature of the Schmidt-Kennicutt relation. 

These results are shown here for the case of a grand- 
design type spiral galaxy but the nature of the shock-driven 
cooling should apply to any galactic driven flow. We ex- 
pect similar processes to occur in flocculent spirals or other 
convergent flows, and hence provide a universal explanation 
for the triggering of star formation and the ensuing rates. 
The necessary condition is of a bi-stable interstellar medium 
where cooling instabilities occur above a critical density. 

The actual values of the star formation rates found 
here are offset to higher values from the extragalactic val- 
ues. Additional forms of support including higher vertical 
turbulence, magnetic fields and feedback could easily lower 
the star formation rates and hence account for the discrep- 
ancy, as would the sudden turning on of self-gravity. Mag- 
netic fields and feedback will delay the on set of star for- 
mation and decrease the s ubsequent rates dBaneriee et all 
120091 ; IPrice fc Bate] 120091 ; IVazquez-Semadeni et al.1 l2010i ). 
while the sudden turning on of self-gravity neglects any pre- 
vious star formation histories in the gas and thus artificially 
increases the rates. Taken together, these effects will neces- 
sitate higher gas surface densities to produce the same star 
formation rates, but given the success here of reproducing 
the Esfr oc Egat, they need only provide a constant offset. 
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